Three-Dimensional Solitary Waves and Vortices in a Discrete Nonlinear Schrodinger 
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In a benchmark dynamical-lattice model in three dimensions, the discrete nonlinear Schrodinger 
equation, we find discrete vortex solitons with various values of the topological charge S. Stability 
regions for the vortices with 5 = 0, 1, 3 are investigated. The S = 2 vortex is unstable, spontaneously 
rearranging into a stable one with 5 = 3. In a two-component extension of the model, we find a novel 
class of stable structures, consisting of vortices in the different components, perpendicularly oriented 
to each other. Self-localized states of the proposed types can be observed experimentally in Bose- 
Einstein condensates trapped in optical lattices, and in photonic crystals built of microresonators. 
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Introduction. Intrinsic localized modes (ILMs), alias 
discrete breathers, in nonlinear dynamical lattices have 
inspired a vast array of theoretical and experimental 
studies. They have attracted attention due to their inher- 
ent ability to concentrate and, potentially, transport the 
vibrational energy in a coherent fashion (see, e.g., [1]). 
Settings in which these objects appear range from arrays 
of nonlinear-optical waveguides [2] and photonic crystals 
[3] to Bose-Einstein condensates (BECs) in optical-lattice 
(OL) potentials [4], and from various systems based on 
nonlinear springs [5] to Josephson-junction ladders [6] 
and dynamical models of the DNA double strand [7] . 

A benchmark model, which generically emerges in the 
description of dynamics in nonlinear lattices, is the dis- 
crete nonlinear Schrodinger (DNLS) equation [8]. It 
finds its most straightforward physical realizations in two 
of the above-mentioned settings, viz., arrays of optical 
waveguides [9,10], and networks of BEG drops trapped 
in OLs [4]. While the former system may be, effectively, 
1- and 2-, but not 3-dimensional, the latter was exper- 
imentally created in three dimensions (3D) as well [11], 
which suggests a direct physical implementation of the 
3D DNLS model. Another physical realization of the 3D 
DNLS equation may be provided by a crystal built of 
microresonators trapping photons [12] or polaritons [13]. 

In spite of the physical relevance of the NLS equation 
in the 3D case, very few works attempted to find localized 
solutions in this system, and those were actually done in 
the absence of OLs [14]. Only very recently, a possibil- 
ity of the existence of stable 3D solitons in continuum 
NLS equations including OL potentials has been demon- 
strated [15]. A problem of fundamental interest, both in 
its own right and in view of the possibility of the exper- 
imental realization (principally in the BEG-OL setting), 
concerns the search for 3D solitons in the DNLS model 
proper. Especially intriguing is a possibility to construct 
stable ILMs with intrinsic vorticity (topological charge), 
which would be an entirely new class of ILMs in 3D. 

In this work, we demonstrate that the discreteness in- 



deed makes it possible to stabilize, in the DNLS model 
with attraction, not only ordinary 3D ILMs, but also vor- 
tex solitons (they are strongly unstable in the continuum 
limit; notice, however, that 3D vortex solitons can be sta- 
bilized in continuum models with competing nonlineari- 
ties [16]). These include not only fundamental discrete 
vortices, with the topological charge S = I, whose sta- 
bility in the 3D case may be surprising by itself, but also 
higher-order vortices, such as ones with S" = 3 (in the 
above-mentioned 3D continuum models with competing 
nonlinearities, stable higher-order vortices have not yet 
been found) . We also extend the considerations to multi- 
component DNLS models, that allow for the existence 
and stability of still more challenging configurations. In 
particular, we introduce a novel type of a compound vor- 
tex, consisting of two vortices with the same value of S 
in the two components, whose orientations are perpen- 
dicular to each other. Such solutions are stable too, in 
certain parametric intervals. 

The Model. The DNLS equation on the cubic lattice 
with a coupling constant C is [8] 
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ILM solutions <f)i^m,n = Gxp{iAt)ui^m,n, whcrc A is the 
frequency (chemical potential in the BEG context) and 
the stationary eigenfunctions ui^rn,n obey the equation 
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Due to the invariance of Eq. (1) against the scaling trans- 
formation, t — > t/U^, C — > CU'^, and u -^ u/U, with 
an arbitrary constant U, one can either fix the coupling 
constant, as C = Cq, and vary A, with the objective to 
explore a full family of solutions of a certain type, or, 
alternatively, fix A = Aq, and follow the variation of C. 
The actual control parameter, that is invariant against 
the scaling, is C/A. 



Solutions to Eq. (2) (generally, complex ones) are ob- 
tained by means of a Newton method. To test their sta- 
bility, perturbed solutions are used in the form [17] 
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where e is an infinitesimal amplitude of the perturbation, 
and A is its (generally, complex) eigenfrequency. The sub- 
stitution of Eq. (3) into Eq. (1) gives rise to linearized 
equations for the perturbation eigenmodes. 
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where Fk = -C{uk+i + Uk-i + Uk+N + "fc-jv + "fc+w^ + 
Uk-N^ — 6ufe)-f Aufe — \uk\ Uk, and the string index, fc, is 
defined so that {l,m,n) ^ k = l + {m-l)N+{n-l)N^. 
Dirichlet boundary conditions were imposed. 

According to the scale invariance discussed above, we 
examine solutions of Eq. (2) by fixing the frequency, 
A = 2, and varying the coupling C (in the BEC context, 
this means fixing the chemical potential, and varying the 
OL strength, as is experimentally feasible). The solutions 
with different values of the topological charge S (here, it 
ranges between and 3) are generated by an iterative 
scheme with an initial ansatz motivated by 3D vortices 
in the continuum limit [16], 
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where (lo,mo,no) is the location of the vortex' center, 
and 77 is a scale parameter. The Newton algorithm was 
then iterated until it converged to 1 part in 10^. Our 
results are typically shown for 9x9x9 and 11 xll xll 
site lattices, but larger ones were also investigated. 
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FIG. 1. The ILM with S = is shown for C = 1. The left 
panel shows the 3D contour plot Ke{ui^m,n) ~ 0.5. The right 
panels show 2D cross sections of the solution through n — 6 
(top) and n = 5 (bottom) for the 11 xll xll DNLS lattice. 

Results. Basic results for the solutions with different 
topological charges can be summarized as follows. Ordi- 
nary ILMs with S = are stable below a critical value 



Ccr ~ 2 of the couphng constant, as they are strongly 
unstable (against collapse) in the continuum limit of 
C —f 00. An example of a stable ordinary soliton is shown 
in Fig. 1. As ILMs with S* = have the largest stability 
interval, C < Ccr , as compared to topologically charged 
ones (see below), they can only be destroyed if unstable, 
rather than being transformed into ILMs of other types. 
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FIG. 2. The top panels show level contours at 
R.e(it;_m,„) — ±0.5 (left) and Im(«;,m_„) — ±0.5 (right) for the 
3D vortex ILM with 5 = 1. Red (color version) / light-gray 
(black-and-white) and blue/dark-gray surfaces pertain to the 
levels with —0.5 and +0.5 values, respectively. Cross sections 
of the vortex are shown in four middle panels, (c) and (d). The 
bottom row displays the development of instability of the vor- 
tex for C = 0.7, through the time evolution of its amplitude, 
and a 2D cut of the profile at i = 100 [(f) top and bottom, 
respectively]. The unstable vortex transforms itself into an 
ordinary ILM with S — 0. The left bottom panel (e) shows 
the spectral plane {Xr,Xi) of the linear stability eigenvalues 
for the same unstable vortex. 

3D vortices with 5=1 have also been found. They 
are stable (see Fig. 2) below a critical value Ccr ~ 0.65 
(similarly to their 2D counterparts [18]). At the instabil- 
ity threshold, a quartet of complex eigenvalues emerges 



from collision of two imaginary eigenvalue pairs (for de- 
tails, see, e.g., Refs. [18,19]). Numerically simulated de- 
velopment of the instability is displayed in Fig. 2, for a 
typical case with C = 0.7 > Ccr • The perturbations 
destroy the vortex structure and, as a result, an ordi- 
nary {S — 0) ILM emerges; obviously, the change of the 
topological charge is possible in the lattice, in which the 
angular momentum is not a dynamical invariant. 
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FIG. 3. The ILM vortex with S = 2 for C = 0.01. Top 
panels have the same meaning as the top panels in Fig. 2. 
The bottom left panel (c) displays the linear stability eigen- 
values, while (d) shows the result of long evolution of this 
unstable vortex. The eventual state, shown through its 2D 
cross-sections at f = 1000, is a vortex with S = 3 (see Fig. 4) , 
which is stable for this value of C. 
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FIG. 4. Stable stationary vortex ILMs with 5 = 3 for 
C — 0.01. Panels have same meaning as top panels in Fig. 2. 

An example of a vortex with topological charge S — 2 
is shown in Fig. 3 for C = 0.01. Similar to its 2D coun- 
terpart [18], this complex solution is unstable through a 
real eigenvalue pair at all values of C (notice, however, 
that purely real stable solutions to the 2D DNLS equa- 
tion, that may be identified as quasi-vortices similar to 
the solitons with 5* = 2, have been very recently found 



[19]; counterparts of such solutions exist in the 3D case 
as well). What is more interesting, however, is that this 
unstable ILM with S = 2 reshapes itself not into the one 
with iS* = or iS* = 1, but rather towards a stable vor- 
tex with S" = 3, as seen in Fig. 3. The stabilization of 
the 3D vortex ILM through spontaneous increase of S is 
a striking result (again, feasible only in dynamical lat- 
tices). In Fig. 4 we show the stable S = 3 discrete vortex 
for the same case, C = 0.01 (stable higher-order vortex 
ILMs were very recently found in the 2D DNLS model 
too [19]). The instability of the vortex with S" = 2 vs. the 
stability of the vortex with 5' = 3 may be understood, 
in loose terms (in the 2D case as well) in terms of the 
lattice-induced Peierls-Nabarro (PN) potential acting on 
the soliton. Indeed, it is the PN potential which may sta- 
bilize a soliton which is otherwise strongly unstable. It is 
seen from Figs. 3 and 4 that the PN potential induced by 
the cubic lattice is, obviously, a much stronger factor for 
the soliton with S — 3 than with S = 2, due to the sym- 
metry difference between the former one and the lattice 
(i.e., the "skeleton" of the vortex lies at angles of 27r/3 
as opposed to the 7r/2 of the underlying cubic lattice) . 
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FIG. 5. A complex of two orthogonal vortices with S = 1 
in the two-component system is shown for C = 0.01. The top 
and bottom panels correspond to the two components, and 
they have the same meaning as the top panels in Fig. 2. 

Another striking feature, unique to the 3D case, is a 
possibility of the existence of vortex complexes in a multi- 
component system, with the vortices in different (up to 
three) components orthogonal to each other. We con- 
sider, in particular, two coupled DNLS equations. 
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which describe an array of BEC droplets composed of a 
mixture of two different species [20] . In the case of the 
model based on the photon or polariton field trapped in 
the lattice of microresonators, and 4' refer to two dif- 
ferent polarizations or distinct cavity modes. In Eqs. (6), 
(3 is the relative, intra-species interaction strength. 

We examine a complex of two orthogonal vortices, in 
which the one in the first component is directed perpen- 
dicular to the {I, m) plane, while in the second compo- 
nent, the vortex is orthogonal to the {l,n) plane. An 
example of such a stable complex is shown, for (3 = 0.5, 
in Fig. 5. We have found that the orthogonal complexes 
may be stable for f3 < 1 (for sufficiently weak coupling C: 
the solutions of Fig. 5 are stable for C < C^. w 0.025), 
and are unstable for /? > 1. This can be qualitatively 
understood in terms of the Hamiltonian of the attractive 
interaction between the two components, each having the 
characteristic "doughnut" [16] vortcx-soliton shape (in 
the continuum limit). Indeed, one can roughly estimate 
the interaction energy (negative) through the volume V 
of the overlap between two cylinders of a radius p (which 
represent long inner holes of the doughnuts) intersect- 
ing at an angle 9, V ^ {20/3)p^/ sin9 (the divergence at 
61 ^ is limited by a finite length of the holes). As it 
follows from here, the interaction energy has a maximum 
at 6* = 7r/2, which corresponds, by itself, to an unstable 
equilibrium state of two orthogonal vortices. Actually, 
the equilibrium is transformed into a stable one by the 
pinning to the PN potential, provided that the interac- 
tion is not too strong, i.e., /3 is not too large. 

Conclusions. The above results are a first step towards 
an understanding of topologically nontrivial ILMs in 3D 
dynamical lattices. Besides generating stable ILMs with- 
out {S = 0) or with {S = 1 and 3) topological charge, 
the 3D lattice gives rise to a variety of novel dynami- 
cal effects, such as the reshaping of the S = 2 unstable 
vortex into a stable one with S' = 3. Furthermore, the 
3D dynamical lattice sustains quite unusual but stable 
states, such as the two-component complex, with the in- 
dividual vortices in the components orthogonal to each 
other. Studying more complex configurations in such a 
rich setup, and examining interactions between ILMs, are 
challenging problems for future work. 

The 3D vortices predicted in this work can be created 
in an self-attractive BEC trapped in an optical lattice 
(OL). Relevant physical parameters are essentially the 
same as those at which this medium is experimentally 
available [4,11,20], and for which various ID and 2D lo- 
calized structures were predicted [4,15], i.e., > 10'* atoms 
trapped in an OL with the period ^ 1 fim and the size 

> 10 X 10 X 10. The soliton-formation time is ~ 1 ms, 
and the experiment can be easily run on the time scale of 

> 1 s. In principle, the same vortex solitons can be also 
created in a cubic lattice composed of microresonators 
trapping photons or polaritons. 
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